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Two-Scale Kirchhoff Theory: Comparison of Experimental Observations With 

Theoretical Prediction 
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We introduce a non-perturbative two scale Kirchhoff theory, in the context of light scattering by 
a rough surface. This is a two scale theory which considers the roughness both in the wavelength 
scale (small scale) and in the scales much larger than the wavelength of the incident light (large 
scale). The theory can precisely explain the small peaks which appear at certain scattering angles. 
These peaks can not be explained by one scale theories. The theory was assessed by calculating 
the light scattering profiles using the Atomic Force Microscope (AFM) images, as well as surface 
profilometer scans of a rough surface, and comparing the results with experiments. The theory is in 
good agreement with the experimental results. 



PACS: 42.25.Fx, 68.37.-d 



I. INTRODUCTION 



Wave scattering by rough surfaces has been extensively 
studied both analytically and experimentally. For analyt- 
ical approaches two methods have been generally consid- 
ered: rigorous electromagnetic theory and approximate 
methods. The Kirchhoff theory is among the electromag- 
netic theories and is known as a "tangent plane theory". 
This theory is most widely used to calculate the distri- 
bution of the specular and diffuse parts of the reflected 
light. The Kirchhoff theory treats any point on a scat- 
tering surface as a part of an infinite plane, parallel to 
the local surface tangent. The theory is therefore exact 
for an infinite, smooth and planar scatterer, but is ap- 
proximate for scatterers that are finite sized, non-planar 
or for rough surfaces [1]. Due to the computational limi- 
tations, most studies have been done for one dimensional 
data of the surfaces. There are only few cases of the 
analysis of two dimensional surface data. One and two 
dimensional exact approaches have been successfully ap- 
plied to dielectric, metallic or perfectly conducting sur- 
faces [2,3], deterministic surfaces [4,5], dielectric films on 
a glass substrate [6] and dielectric films [7,8]. Such ex- 
act calculations have been compared with experimental 
results and approximate models [6,9]. Also some authors 
studied wave scattering from random layers with rough 
interfaces [10,11]. 

The joint probability density functions (PDF) of sur- 
face slopes and heights P(d x h,h), is a key function in 
the estimation of the main parameters of wave scattering 
by a rough surface [12-16]. This is more obvious in a 
geometrical optics approach, when the angular distribu- 
tion of the scattered power is proportional to the specular 
reflecting slope PDF. The slope PDF has also been in- 
troduced in references [12,15,16] in the context of Bragg 



scattering. They have shown that, the Bragg scattering 
results must be averaged by the proper slope PDF of the 
rough surface. This is also true for the estimation of the 
thermal emission from rough surfaces at small grazing 
angles [15,16]. 

In the present paper, we introduce a non-perturbative 
two scale Kirchhoff theory. The theory is applied to ex- 
plain the small peaks observed in the scattering profile 
of a rough surface, at certain scattering angles. The the- 
ory employs the data obtained from the rough surface in 
two different scales. To check the theory we have mea- 
sured the scattered light intensity MS M. function of the 
scattering angle, 1(0), using a setup consisting of a He- 
Ne laser (632. 8nm), a photo-multiplier tube (PMT) de- 
tector and a computer controlled micro-stepper rotation 
stage. The resolution of the micro-stepper was 0.5 min- 
utes. Alumina sheets were used as the rough samples. 
The surface topography of the alumina samples in small 
scale (< 5/im) was obtained using an atomic force micro- 
scope (AFM) (Park Scientific Instruments). The images 
in small scale were collected in a constant force mode 
and digitized into 256 x 256 pixels. A commercial stan- 
dard pyramidal S^A^ tip was used. A variety of scans, 
each with size L, where recorded at random locations on 
the surface. The large scale (< 5mm) morphology line 
scans of the alumina samples were recorded using a sur- 
face profilometer (Taylor Hobson). Figures (1) and (2) 
show typical AFM image and surface profile data with 
resolutions of about 20nm and 0.25/im, respectively. 




FIG. 1. AFM image of the Alumina surface in the length 
scale 5/xm x (small scale). 

II. NON-PERTURBATIVE TWO SCALE 
KIRCHHOFF THEORY 

The Kirchhoff theory is based on three major assump- 
tions [1]: 

a) - The surface is observed from far field. 

b) - The surface is regarded as fiat, and the optical be- 
havior is locally identical at any given point on the sur- 
face. Therefore the Fresnel laws can be locally applied. 

c) - The amplitude of the reflection coefficient, i?o is 
independent of the position on the rough surface. 

The field scattered by the rough surface, ip sc (r), is ob- 
tained by an integration over the mean reference plane 
Sm [1], (the geometry is displayed in figure (3)) 

ike^kr) I I { M +b M_ c) 
4?rr J J SM dx dy 

exp (ik(Ax + By + Ch(x ,y )))dx dy (1) 
where 

A = sin 6*i — sin 6*2 cos 9%, 
B = — sin 82 sin 6*3, 
C = —(cos 6*i + cos 62), 

a = sin6>i(l — i? ) + sin 2 cos 8 3 (1 + R ), 

b = sin 6» 2 sin 6 3 (1 + R ), 

c = cos6> 2 (l + Ro) - cos6»i(l - i? ) 

In the derivation of the equation(l), it is assumed that 
the incident wave ip m is a plane wave with a wave vector 
k as ip m (r) = exp(ik • r). 
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FIG. 2. Profilometer scans of the Alumina surface with res- 
olution 0.25pm (large scale). 

In most cases, the wave scattering models from rough 
surfaces implicitly assume that the surface is rough on a 
single scale. However, in practice all surfaces are rough 
on several scales, ranging from atomic scale to the scale 
determined by the length of the surface. Nevertheless, 
only a finite range of scales are important in scattering 
of waves from a surface, i.e. the range covering the wave- 
length of the incident radiation. Models have been devel- 
oped for describing surfaces that consist of high frequency 
fluctuations superimposed on a slowly varying roughness 
[1]. These models use perturbation theories to describe 
the scattering from the high frequency roughness and this 
is modified in some manner by the low frequency com- 
ponent [17]. All of the perturbative methods deal with 
the effect of the large scale fluctuations as perturbation to 
the small scale height fluctuations. Here, we intend to ob- 
serve the surface in two scales with resolutions of nanome- 
ter and micrometer. The figure (4), shows schematically 
the modulation of small scale height fluctuations by large 
scale variations. Various statistical parameters like the 
joint height and height gradient PDF, surface roughness 
<7, correlation function C(i?), correlation length r etc., 
were measured in two scale. 

In what follows, we are going to describe the non- 
perturbative two scale Kirchhoff theory. We first cal- 
culate the contributions of the coherent and the diffuse 
fields by the Kirchhoff theory in small scale. The coher- 
ent field with a gaussian height distribution will be [1] : 

<^ c ><^ sc >*=/ cxp(-<7) (2) 

where g — k 2 a 2 C 2 . Also fc, a and Iq are the norm of 
wave vector, surface roughness in small scale and the 
scattered reflected intensity of the corresponding smooth 
surface. For isotropic surface and for samples with the 
sizes much larger than the correlation length L»t, (and 
for a slightly rough surface i.e. g <§; 1) the diffuse field 




FIG. 4. Two scale observation model of Alumina surface. 

intensity for Gaussian height distribution will be given 
by [1]: 

< Id >= 

fc 2 F 2 r 2 . . A . k 2 (A 2 + B 2 )r\ 
47rr 2 9 cxp(-ff) A M cxp( ) (3) 

where, F = ^(4r + "fr + an< ^ ^ s the effective area 
of rough surface which experience the incident radiation. 
Therefore, the overall scattered intensity is written as [1] : 

< I >= I cxp(- 9 ) + <I d > (4) 

So far, we have expressed the results of the light scat- 
tering from the surface in small scale. Now we divide the 
whole surface to many small pieces (meshes) such that 
length of which is smallest scale of our observation. In 
each mesh we can apply one scale ( small scale) Kirchhoff 
theory. Therefore for each mesh we have a similar expres- 
sion for coherent field as equation (2), but with different 
angles which depend on the positions of the small mesh. 
In the small scale, we denote the height field in position 
x s and y s with h s . Therefore, one can write the height 
field in any position, x, as follows: 

h= hi + h' s 
x= xi + x' s 



FIG. 5. A two dimensional scheme, shown the necessary 
parameters, h' s , here is the h s rotated such that it coincide on 
the vertical time In the three dimensional case, h' s is rotated 
to h a by Euler matrices. 



y= vi + v's (5) 

The indices I and s denote the large and small scales, 
respectively. The vector ( x' a , y' s ) is the position of 
the h' s on the small scale coordinates. In figure (5) 
we have shown the {h! a , x' s , y' s ) and (h s ,x s ,y s ), schemat- 
ically. We note that the AFM images will gives us the 
h s (x s , y s ) and via the large scale topography we will find 
the hi(xi,yi). The vectors (h' s ,x' s ,y' s ) and (h 8 ,x s ,y s ) can 
be related to each other, via rotational Euler matrix, 
with three rotational angel a, (3 and 7 i.e. A(a, /3, 7) = 
R h (~f)R y {f3)R x (a). 

The Local angles (#1, O2, #3) are defined by the average 
plane in the small scale. Therefore, all a, b, c, A, B, C are 
constant for all points within the small piece. In each 
small scale element hi is fixed so that J-j^ 1 - = 0. Hence, 
the total scattered field has the following expression: 

-c s )expik(A s x' s + B s y' s + C s (h' s {x' s ,y' s )) 
dx' s dy' s ] expik(A t xi + B t yi + Qhi) 

= ^ SC M eM*k(Aixi + B m + Qht)) (6) 

%i,yi 

we note that = -S^- and the summation is over the 

OX B OX s 

small scale samples modulated by the large scale fluctu- 
ations. We assume that the joint PDF of heights and its 



slope of two scales are independent, then the average of 
the field scattered in any direction will be given by: 

< v-Ar) >=n y E <^ c w > 

hifixhi xi ,yi 

exp(ikCihi) cxp(ik(Aixi + Biyi))P(hi,d x hi) 
= N Y <^ c (r)>cxp(ik(C l h l ))P(h l ,d x h l ) 

hi,dxhi 

Y exp((ik(AiXi + Biyi))) 



N 



sm(kL x ) sin(kL y ) 



kLrr 



kLii 



A 



M 



Y <r s c (r)>cxp(ikC l h l )P(h h d x h l ) (7) 

hi,d x hi 

The subscript (— e) denotes scattering from the sur- 
face without the edge terms. The rough surface has 
been assumed to be rectangular with extent —X < xo < 
X, — Y < yo < Y. Also L x and L y are length scales in 
the scattering area (the effective area of light incidence), 



and S a 



siYi(kL x ) sin(hLy) 



Am is the constant term 



'A/ - — kT x kL^^M 

all observation angles. The quantity NP(hi,d x hi) is the 
number of points with height h and slope d x hi . It is noted 
that for a homogeneous surface p(h, x) is independent of 
position along the surface, x. In order to do analytical 
calculation, it is necessary to assume that the edge effects 
are non-stochastic, i.e. < ip e >= V J e [1]- Based on this 
assumption, the coherent part becomes: 

< i coh >=< r c >< r c >*= 

N 2 Y I < ^ > \ 2p (hi,d x hi)P(h 2 ,d x h 2 ) 

hi ,d x hi Ji2,dx h-2 



exp(ik(C(h 2 - h^)) 



(8) 



where ip sc = ip e + ip- e . It is noted the non-stochastic 
assumption of the edge effect leads to the cancelation of 
all terms containing edge effects. In cylindrical coordi- 
nates, for an isotropic surface, the substitutions X2 — X1 = 
RcosO and y 2 — Vi = RsinO can be made. Since the 
heights PDF and the heights difference PDF are inde- 
pendent (we will confirm this assumption in the next 
section), i.e. P(h\,d x hi) = P(h\)P(d x h\). 

Define, 

Y dhidh 2 exp(ik(C(h 2 - hi)) P (hi) P (fa) 

hi,h2 

= X (kC,~kC, R) 
then one finds: 

< i coh > = s 2 M \ y xp(d x hi) < r s c > 1 2 




P(h)p(d x h) 

FIG. 6. Joint PDF vs P(h)P(d x h), that shows the height 
and slope PDFs are almost independent. 



x(kC,-kC, R) 



(9) 



It is known that the total average scattered field in small 
scale is < ipf >= x{kC s )ip s c . 

For a gaussian height distribution, the one and two- 
dimensional characteristic function is given by: 



x(kC s ) 



1 f +co h 2 
7== / cxp(-— — )cxp(ikC s h s )dh s 



exp(-^C> s 72), 
and 



(10) 



X (kC, -kC, R) = cxp(~k 2 C 2 a? (1 - C(R))). (11) 

where C(R) = <h ^ r ) h ^+ R )> ; j s the surface correlation 

function in the large scale. Also the average of total 
intensity are given by : 

< i tot >=< r c r c * > 



N E E E p (hi,d x hi) < r s c r/ > 

x ,yo xi,yi 



exp(ik(A(x 2 - x x ) + B(y 2 - yi))) 



(12) 



Performing the summation we find Y] h — N, where N 
is the number of points on the surface. So, the average 
total intensity becomes: 

< i to t >= s 2 m n y N(o x hi) < r s c r s c " > (13) 

Finally, the diffuse field intensity is obtained as: 

< I d >=< Rot > - < Icoh > (14) 
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FIG. 7. Comparison of theoretical prediction via two scale 
Kirshhoff theory and experimental results for scattered field 
(bold symbols). 

III. COMPARISON WITH EXPERIMENTS 

Here we test the non-perturbative two-scale Kirshhoff 
theory with experiment. For this purpose, we obtain 
the height profile of the Alumina sheets as the rough 
samples, using the profilometer in large scale and the 
AFM images in small scales. Indeed we intend to ob- 
serve the surface in two scales; nano-meter and micron. 
To use the two scale theory the surface must possess 
two conditions. First, the PDF of the height and its 
slope must be independent at small and large scales, 
i.e. P(hi,d x h h h s ,d x h s ) = P(hi,d x hi)P(h s ,d x h s ). The 
homogeneous rough surfaces possess this condition. In- 
deed statistical parameters in small scale ( roughness, 
exponents, etc.) are similar at any point of the sam- 
ple (large scale). This means that the two PDFs are 
independent. The second condition is that the height 
and height gradient fluctuation must be independent 
in the large scale. This means that the joint PDF of 
the height and height gradients can be decomposed as 
P(hi,d x hi) = P(hi)P(d x hi). This assumption needs con- 
firmation. In figure (6), we have plotted the joint PDF 
P(hi,d x hi) vs P(hi)P(d x hi). It is obvious that the joint 
PDF vs multiplication of single PDFs fits with a line 
with slope one. Considering its statistical error we ob- 
serve that the height and height gradient PDFs are inde- 
pendent. For large values of h and d x h, our assumption 
becomes poor and thus uncertainty increases. 

To compare the experimental observation with those of 
the theoretical prediction, we need estimate the several 
types of PDFs in small and large scales. In the equa- 
tion (8), we need to evaluate the quantity < ^f c -0f c > 
in the small scale and PDF of the height gradients in 
the large scale. To evaluate the intensity < tp s s c ip s s c >, 
we have to use the equation (2), where the averaging is 



done in the small scale. Therefore we need the PDF of 
height fluctuation in the small scale. Also we need other 
statistical quantities such as surface roughness cr, corre- 
lation function C(i?), correlation length r etc., in small 
and large scales. We evaluate the height -height corre- 
lation function < h(x + R)h(x) > vs radial distance R 
for large scale fluctuations. We find the following expres- 
sions for the Alumina surface as, C(R) = 2.14 exp(— 
and 1.27 exp(— 0.58i?), for small and large scales, respec- 
tively. Also the roughness exponent, variance and scal- 
ing length for the small (large) scale have been found as, 
0.85 (0.85,031), 0.31/mti (1.33^m) and 1.5^m (19.4/jm), 
respectively. It is found that the height PDF in the 
two scales are almost gaussian. The estimated statistical 
quantities enable us to predict the average total inten- 
sity. In figure (7), we have plotted the experimental ob- 
servation and theoretical prediction of total intensity. It 
is evident that the theoretical prediction fits with those 
of experimental observation. We observe the theory is 
able to predict small peak in the angle ~ 18° in the 
variation of the total intensity vs angle scale 02- We 
note that if one plots the PDF of height gradient, then 
finds that the PDF has also small peaks at angle scale 
tan -1 ^/).) ~ 9°. This means that the gradient PDF is 
responsible to have a small peak in the variation of the 
total intensity in terms of angle scale (we note that the 
slope a = tan~ 1 (d x h), we produces 2a contribution in 
the reflection of the light from the surface). In figure (8), 
the behavior of the slope PDF (d Xl hi) in terms of d xi hi, 
has been given. Also as shown in figure (7), the two scale 
Kirshhoff theory is able to predict the small peak in the 
variation of the total intensity in terms of angle scale. 
As we observe, there are other peaks in the figure (7), 
where the theory can not predict the peaks for large an- 
gle scales. Indeed for these angle scales we should take 
into account the shadowing effect [15,16]. In ref. [12], the 
validity range of geometrical shadow functions has been 
investigated for a randomly rough surface for which the 
shadowed Kirchhoff approximation has been shown to 
give good results for the scattered intensity distribution. 
We will discuss the modification of the two scale Kirsh- 
hoff theory by the shadowing effect elsewhere. 
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